Laser Plasma Interaction and Non-classical Properties of Radiation Field 
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We show by exphcit calculations that non-classical states of the radiation field can be produced by 
allowing short term interaction between a coherent state of the radiation field with plasma. Whereas, 
long term interaction, which thermalizes the radiation field, can produce non-classical states of the 
radiation field only at sufficiently small temperatures. A measure of fc*'^ order squeezing, stricter 
than the one proposed by Zhang et al [1], is used to check the emergence of squeezing. It is also 
shown that photons in the considered thermalized field would follow super-Poissonian statistics. 
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I. INTRODUCTION 

Most of the light sources around us are classical in 
nature; in the sense that all the properties of the emit- 
ted radiation can be explained by assuming that the 
emission process in the source is actually some clas- 
sical stochastic process. To put in more pedagogical 
terms, a radiation field is said to be 'classical' if the 
Glauber-Sudarshan P-function [2] [3] [4] corresponding 
to the state of the radiation field, is positive definite. 

It turns out that there can be states of the radiation 
field for which the 7^-function is either highly singular 
(for example, containing derivatives of the Dirac delta 
function) or is not positive in the whole parameter space 
(sometimes also referred to as the phase space). Prop- 
erties of such radiation fields cannot be explained by 
modeling the emission process by any classical stochas- 
tic process. It is with these light sources that one can 
exploit the full power of the quantum nature of the elec- 
tromagnetic field. Some of the properties which can 
only be exhibited by 'non-classical' radiation are an- 
tibunching, sub-Poissonian distribution of the number 
of photons in the field (sub-Poissonian statistics) and 
squeezing. 

One way to produce non-classical radiation is via the 
interaction of a coherent radiation (easily obtained from 
a laser operating well over threshold) with free electrons 
[5]. In this work we calculate the short term and long 
term[6] non-classical effects of interaction, on the incom- 
ing coherent radiation, with the free electron gas. We 
explicitly show that short term interaction can be used 
to produce non-classical states of radiation field exhibit- 
ing squeezing. And contrary to expectation, even long 
term interaction, which arguably decoheres the radia- 
tion field, can be used to produce non-classical radiation 
at sufficiently small temperatures. 



II. FIELD HAMILTONIAN IN PRESENCE OF 
PLASMA 

Ben-Aryeh and Mann [5] have shown that under cer- 
tain approximations the Hamiltonian for the pair of 
modes, k, A and — k, A of the radiation field interact- 
ing with slow moving and sparse 'free' electrons (which, 
from here on, we will refer to as plasma) in a volume 
having dimensions larger than |k| can be written as 



(1) 



where 



a- = (ak,A - a-k,A) /V2. 



(3) 
(4) 



is obtained from by replacing all the -I- in 
the subscript by — . uj and f2i are real constants re- 
lated to ujk — ck and the density of free electrons in the 



plasma, and il2 is a complex number. 



and a_ are 



bosonic annihilation operators for two uncoupled modes 
of standing waves. 

For our calculations we start with a more general 
mode pair Hamiltonian [7], 
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n2a^ + nia (5) 



cli is a real constant, whereas, fli and fl2 are con- 
stant complex numbers. One can easily identify 
1/2 (a'l'a -I- 1/2) , <^^^/2 and °-^/2 in eqn(5) with the elements 
of the SU(1,1) algebra, K3, and K^, respectively 
(harmonic oscillator realization of SU(1,1) algebra). K^, 
K+, a and together form the elements of the 
double photon algebra [8]. Hence, 



H = 0JK3 + OiA'+ + rj*A'_ 



2a' 



(6) 
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This is the most general Hermitian operator that can 
be constructed from the elements of the double photon 
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algebra, modulo a real constant. Eqn(6) can be rewrit- 
ten as 



where 



a 

a = I t 



"/2 ni 



(7) 



(8) 



(9) 



(10) 



H can be simplified to an element of the SU(1,1) al- 
gebra (apart from a c-number) by the unitary transfor- 
mation described below. For 



a ■ 



(11) 



H' = D{a)HD^{a) =P-c (12) 
where, D(a) is the displacement operator defined as 

D{a) = exp (aa^ — a* a) (13) 



A. The Disentanglement Formula 

To deal with coherent states it is best to write the 
SU(1,1) group element in eqn(16) in the following dis- 
entangled form, 



(17) 



That we can always do this, ie, such P, 7 and 6 would 
always exist is guaranteed by the Wei-Norman decou- 
pling theorem [9]. To find 7 and S in terms of co and 
fli we consider the following relations. 



-it{ad P)q^ 



(18) 



where the super operator {ad P) is defined recursively 



as 



and 



{ad PTA = [P, A]n = [P, [P, A]n-i] (19) 



P = 



One can show that 



-ttp _ / cos^t + i^sin^t 
e - I nt 



j% sin^t 



i-r- sin (l)t cos 0t — i-^ sin (j)tj 



(20) 



(21) 



P = -aUa = UK3 + aiK+ + niK_ (14) 

and 

c = OU-^0. (15) 



III. UNITARY TIME EVOLUTION 

The state of the radiation field would undergo unitary 
time evolution for interaction time small compared to 
the dccohcrcncc or relaxation time scale of the system. 
For such small interaction times the unitary time evolu- 
tion operator for the radiation field is given by (7l = 1) 

U{t) = exp{-iHt) 

= D^i {a) exp {-iH't) D(a) (16) 
= e*'=*Dt(a) exp(-iPf)D (a) . 

The initial state of the system under consideration is a 
coherent state. If fii = 0, a coherent state would remain 
coherent at all times and no non-classical properties can 
be observed. For our system f^i ^ 0, but owing to the 
small density of free electrons in the plasma, <C 
w. Nevertheless, in the sections to follow we will show 
that even a small fii leads to evolution of non-classical 
properties in the radiation field. 



where (j> = ^'^^/^ — Using eqn(21) in eqn(18) and 

comparing it with 



yay-i = e""/^(a-^at) 



7 = -2 In 



UJ 



cos^t -|- i— sin (j)t 
2(p 



iVti sin (j)t 
) cos (/it -I- i% s,m(j)t 



zf2J sin (j)t 



)COS(j)t + sinc^t 



(22) 



Faty-i = (eV2 + Se"'^')a^ + Se'^^'a (23) 
one obtains 



(24) 



(25) 



(26) 



Eqn(24), (25) and (26) imply that apart from the physi- 
cally irrelevant phase factor (e*''*) the evolution operator 
for the radiation field is periodic in time, with period, 
T = 27^/0. If r is small compared to the decoherence time 
scale of the system, this periodicity will be refiected in 
all the physical properties of the radiation field. 
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B. Non-classical Radiation 

As discussed in the Introduction, a state of the radia- 
tion field is said to be 'non-classical' if the corresponding 
Glauber-STidarshan "P-function is either negative over 
some region of the parameter space or is more singular 
than the Dirac-delta function. With eqn(16)-(17) and 
(24)- (26) in hand we can calculate the 'P-function for 
the time evolved coherent state to check whether it can 
show non-classical properties or not. 

Let initially the field be in the coherent state |A). 
Therefore, 



p{t) = U{t)\X){X\uHt). 



(27) 



We also have, 



pit) = I d^rj\v){v\V{v,t) (28) 

as the defining relation of the 'P-function. This implies, 
{-mm = Id'v'P{v,t){-C\ri){v\C) 

(29) 



or. 



= d\{-C\p{tme\^^'ei'^-^^'. (30) 

Using eqn(16),(17) and (24)-(26), one can easily eval- 
uate 



ixium) 



g7/4gictg-i(|xr + ICl') 

X e [5 (;8x*'+'5C'+2e''/^x* C+PX* +q<:+r)] 



where 



(31) 



p = 2 {e'/-' - 1) a + 2l3a* 

q = 2 {e'/' -l)a* + 25a (32) 

r = pa*"^ +5a'^ + 2{e'l^ -l)\a\^. 



Therefore, 

(-CIpWIC) 



{-C\u{t)\\){x\u{t)\c) 

|eV4|%-(lCP+|Ap)g[i(/3C*'+/5*C^)] 

^^yi{{2e-"-X+p)C- {2e^'^X+p) * C}] 
^g[iie(5A"+gA+r)]_ 



(33) 



Using eqn(33) in eqn(30) 

V{ii, t)=cj d2^eK/3C*^+/3-C^)eVC*-V*C (34) 



where 

C = 
v' = 



V 



^ JeVf e-l^l^ exp [Re {6X^ +qX + r)] 



Writing <^ = (1 + <2, ;0 = /3i + 1^2 and rj' = r][ + 
where Cfcj/^fe and r?^ (fc = 1,2) are real, one obtains 

V{v,t) = C J dCidC2e(^^^'-''^^^'+'''=^^^^+''^^''^-''f^''0 

(36) 

which blows up for all values of r] (because /3i is either 
> or < 0). Note that ^ is a periodic function of time. 
Also, f3{t = 0) = implies that the 'P-function evaluates 
to 5^(77— A), as expected for the initial coherent state |A). 
This check ascertains the correctness of our calculations. 

On the basis of above calculations we can safely con- 
clude that on brief interaction with plasma a coherent 
state of the radiation field can evolve into a non-classical 
state. In what follows, we explicitly demonstrate that 
the field does show the non-classical property of squeez- 
ing. 



C. Squeezing 

The measure of first order squeezing, as defined in the 
Appendix A, 



where 



I?i = i((ata)-|(a)p-|(Aa)^|) 



{Aaf = {a') - {af. 



Now, 



(ata)-|(a)|2 = {X\W {t)a^ aU it)\X) 

- \ {X\Ul{t)aU{t)\X)\ 



(37) 
(38) 

(39) 



and 



{AaY 



where 



{X\W{t)a^U{t)\X) - {{X\W{t)aU{t)\X)y 
f\\9\- 



f 



cos (j^t + sin ( 



g = ^sin0t (cfeqn(21)). 
Using eqn(39) and (40) in eqn(37) 

'D, = l\9\{\g\-\f\). 

As mentioned earlier, <^uj. Therefore, 
Vi < 0. 



(40) 
(41) 

(42) 
(43) 



This shows that on short term interaction with plasma 
the radiation field becomes 'non-classical' and shows 
first order squeezing. Note that Viit = 0) = 0, as ex- 
pected for a coherent state of the radiation field. Also, 
if the decohcrence time of the system is much larger 
compared to '^'^/<p (period of d and e) the radiation field 
would return to the state of zero squeezing periodically. 
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IV. THE THERMALIZED RADIATION FIELD 

On long term interaction (interaction time greater 
than the decoherence time scale of the system) the field 
comes to thermal equilibrium with the plasma. Con- 
sidering the entire system to be in contact with a heat 
bath at constant temperature; T, the; time independent 
density matrix for the radiation field is given by 



where 9 = ksT and 



Z = Tr (e-"/^) 



(44) 



(45) 



We now exploit the connection between the unnor- 

malized density matrix (e""/") and the time evolution 
operator U(t). The former is obtainable from the later 
by analytically continuing the time parameter to the 
purely imaginary value of —i9~^. Thus from eqn(16) 
the unnormalized density matrix 



-"^0 =e^'D^{a)e-^^'D{a). 



(46) 



Also, following section III A, we can immediately write 
down the disentangled form for the SU(1,1) operator in 
eqn(46), where the disentanglement coefficients are now 
given (by replacing t with —i9~^) by 



one finds that (— C|p|C)e'^'^ is a Gaussian with the coef- 
ficient matrix for the quadratic part being 



M 



V2 _ ;3/2 _ 5/2 -/to(/3) 



-/m(/3) 
Clearly, Tr{M) > 0, and 



Det{M) = eT-|/3|2 

= e^(l-^sinh^l) 



(51) 



(52) 



Of course, if ili = 0, this is always positive. However, 

though \ is very small (compared to cj) in systems of 
physical interest, it is easy to see that for temperatures 
smaller than 



oj/ [21n(a;/0i)] = 



(53) 



Det{M) will become negative. At these temperatures, 
the Gaussian form will have a growing exponent, and 
thus the ■p-function would cease to exist. Thus, in con- 
trast to the thermalized free field (cf fli = case), we 
see that the thermalized radiation field in laser plasma 
interaction can become non-classical at sufficiently low 
temperatures. In the following subsection we explicitly 
show that the thermalized radiation field in this case 
does exhibit the non-classical property of squeezing in 
the above predicted temperature range. 



7 = -21n 



/? = - 



cosh — -I- — sinh ■ 
9 20 



fli sinh i 



' cosh ■ 



+ I sinh I 



ni sinh I 



) cosh ■ 



+ % sinh ■ 



/3*. 



(47) 



(48) 



(49) 



It is now a matter of algebra to check the existence 
of the 7^-function. From eqn(30), we know that the 
■P-function is the two dimensional Fourier transform of 

C|p|C)e'''' • If we can show that this is a well behaved 
function of C,, the existence of the corresponding V- 
function would be ascertained. If, however, (— C|p|C)e'^' 
is not well behaved for some range of 9 we can expect 
to observe non-classical properties of the radiation field 
in that temperature range. 

Z{-C\p\C) = e'o{-Q\D^a)e-'/oD{a)\0 

= eV«+V4e[-H2(i+e"'/')lCI^-^C*^-5Cn] 

4{a(l-e-/^)-/3a*K*-{a*(l-e^/^)-5a}C] 



■i{2(l-e^/^)|c,|^-/3a'' = -5a=}] 



where we used the disentangled form of e ^i" . Now, 
writing as (^i + iC,2-, where Ci and C2 are real variables. 



A. Squeezing 

We will use the disentanglement relation first to 
calculate the partition function in the presence of 
'sources'. That is, we will replace the e^^/" in cqn(46) 

by e^("^-"*)e-''/''e''(''-"). Note that we have not in- 
troduced the source terms as addenda to the Hamilto- 
nian, but rather have put them in as disentangled terms. 
This has the advantage of simplifying the relationship 
between the partition function and various mean val- 
ues that we need to calculate to evaluate Vi and X'2 
eqn(A8). Thus, with the unnormalized density matrix 
with sources denoted by p{s,r]), we find 



Z{e,rj) 



Tr 
Tr 



e°/'Z)t(a)e^(''*-"*)e- 



This implies 



P/egr,(a-a)£,(-Q,) _ 

(54) 



1 am+n 

(aV") = ^ Z(e,r]) 



£=77=0- 



(55) 



Shifting the relative position of the source operators 
with respect to the factor e"^/" will give us a similar 
formula for the mean values of normally ordered op- 
erator products. However, we will soon see that the 
advantages we gain in using the coherent state basis in 
calculating the partition function makes the form used 
most convenient. 
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Using the coherent state basis 



Z{e,r]) = ^ y"d2^(C|e^(«^-«*)e-^/^e''(«-«)|C). (56) 



Clearly, Vi goes negative for 9 < 6c (eqn(53)). This 
is in agreement with the temperature range where the 
■P-function is not defined. 

We note here, that if we follow Zhang et al's [1] defi- 
nition of squeezing parameter 



Where we used the cyclic property of TrQ. Now, using 
the disentangled form of e~^/^ 



V8 + T-/4 



X exp 



2{(l-eV^)-05} 



— sa — rja 



(57) 

This is a simple Gaussian in 77 and e. This allows us to 
calculate the various expectation values simply in terms 
of the parameters /3, 7 and 5. Denoting the deriva- 
tives of the logarithm of the partition function by 

(for (^).=,=o)'*-. (for (^)^ 
find that the only non-zero fs are 



' £=J7=0 



etc., we 



f* — _ ill „„f}, A. 



t^r, = ^(l + ^COth^) 



(58) 



One can now easily show that the following relations 
hold 



(a) 
(aat) 



Tin (-^ 



vn 



/4 _|_ , 0+2 



(59) 



Using eqn(37) 



2^1 = lit 



2 y^'^v 
1 



\A4, 20 ) 



1) 

coth ^ - 



(60) 



Zhang 



2(ata) + (a2) + (at^) 
ter, + Re{tnn) - 1] . 



(62) 

For the physical situation analysed by Ben-Arych and 
Mann [5], Qi being a real positive quantity and Qi <^ co, 



Zhang 



UJ (f> 1 

4^ -2^]'=°*^ 2^ "2 



(63) 



Which would predict squeezing in the same temperature 

range as predicted by the measure of squeezing used by 
us. But, for situations where i?e(Oi) < 0, 



Zhang 



1 

> - 

- 4 



^°^^2^ 



1 



> 0. 



(64) 



That is, Zhang et al's definition of the first order squeez- 
ing parameter does not predict squeezing in such situ- 
ations. Hence, in general the definition of first order 
squeezing used by us is more sensitive than the one pro- 
posed by Zhang et al. In Appendix A (eqn(A9)) we 
argue that the k^^ order squeezing parameter defined 
there (which we have used in our calculations) is a better 
measure of A:*'^ order squeezing than the one proposed 
by Zhang et al; 



Zhang 



(65) 



Using eqn(65), one can make an useful comment 
about the statistics of photon number in the thermalized 
radiation field. Firstly we note, 

(An)'-(n) = \t,,\^ + {t,^-lf 

-F2|a|2 {t,^ -1) + [tesa'^ + t^r^a*^) 

2])Zhang 

(66) 

where n = a^a and fli is taken to be real to arrive at 
the last equality. Next, 



Since, for physical situations l^il <C w. 



2|r!i| 



coth ■ 



49 



(61) 



[i^^ + 2\a\Hsr, - 2\a\^ + 1- 2ter, - \2aHnn + i: 



Note that tnr, is real. Therefore, 



2 

7777 1. 



(67) 



[t% + 2\aWr, - 2\a\^ 



T>2 > 



(j^er) ~H ^7777) i^er) ^-q-q) ^1*^1 if'S'q ^7)7]) ^\^\ ^^er) H~ l] • 



(68) 



Revoking Oi <C w 



because > 1- Prom eqn(65), (66) and (69) 



^2 > 

[tl^ + 2\aWr, - 2|a|2 - 2t,r,j- 1] = (eg) 



{Any - {n) > 0. 



(70) 



\t,r, - if + 2|a|2 {t,r, - 1) 



> 
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Which implies that the photon number in the thermal- 
ized radiation field follows super-Poissonian statistics; a 
property usually exhibited by 'classical' radiation fields. 



oral results, including the definition of the new measure 
of k^^ order squeezing described and used above, and 
some text have been borrowed from [10] with ADG's 
permission. 



V. CONCLUSION 



On the basis of this work we arrive at the following 
important conclusions. Firstly, we found via explicit cal- 
culations that unitary time evolution of the state of the 
radiation field under the Hamiltonian in eqn(5) would 
lead to the production of non-classical states having in- 
definite P-functions. Among the observable properties, 
these non-classical field states would exhibit first order 
squeezing (as also found by Ben-Arych and Mann [5]), 
with the amount of squeezing depending periodically on 
the interaction time between the field and the plasma. 

Secondly, we could also show that even the ther- 
malized radiation field state will have an indefinite V- 
function below a certain temperature {6c). Also, using 
a new measiirc of fc"^ squeezing, expounded up on in 
Appendix A, the thermalized radiation field is shown 
to exhibit first order squeezing at temperatures below 
Oc- The authors feel that it would be interesting to take 
a closer look at exactly what is happening at temper- 
atures close to 0c, so that, contrary to the belief, even 
thermalized field starts exhibiting non-classical proper- 
ties. 

Here we must point out that for the case Re{fli) > 
in eqn(5), the measure of k*^ order squeezing proposed 
by Zhang et al [1] does also predict first order squeez- 
ing in the considered thermalized field, for temperatures 
below 9c- But, it does not predict squeezing in either 
first or second order for i?e(Oi) < case. 

Finally, using cqn(A9), we establish that the con- 
sidered thermalized radiation field would follow super- 
Poissonian statistics at all temperatures. 



Appendix A: fc* order amplitude squeezing 

We define fc"^ order quadrature operators as follows, 

X,{9) = ^{a''e-'' + a^''e''). (Al) 

The k*'^ order quadrature operators used by Zhang et al 

[1] to define k^^ order squeezing are Xfe(O) and Xj.('^/2). 
More generally, one can use Xk{9) and Xk{0 + '^ji) to 
define A;*^ order squeezing. 

[Xfe(^),Xfe(^ + V2)] = |K'at'=] (A2) 

which implies 

(AXfe(0))2(AXfe(0 + V2))' > ^\{{a\a^^\)\\ (A3) 

The quadrature X]^{9) is said to exhibit fc"^ order 
squeezing if 

(AX,(e))^<iK[a^at'=])|. (A4) 

Let 
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1 [(at'=a'=)- |(a'=)|2 + |Aa'=|2cos2(0-0)] if [a^at'=] >G 
Vk{e) = { (A6) 
i [(a'^at'^) - |(a'=)|2 + jAa'^p cos2(<^ - 9)] if [a^at'=] < 

where cj) = Ar5((Aa*). Note, Zhang et al [1] did not consider the case of [a*^,a^*^] < 0. Incorporating this case in 
their definition of squeezing parameter, 

r i[(at'=a'=)-|(afc)|2 + |Aa'=|2cos2(/>] if[a^at'=]>0 

^jZhang = = 0) = (A7) 

\ i [(a'=otfc)-|(a*=)|2 + |Aa*=|2cos2(^] if [a'=, a^^] < 0. 

According to Zhang et al [1], a field exhibits fc**^ order squeezing if p^^^^s < q. 

Using the liberty of choosing 6 in eqn(A4) we arrive at a stricter definition of squeezing. We will say that a field 



is squeezed to order kiiVk < 0, where 

r i [(at'=a'=) - |(a'=)|2 - jAa'^p] if [a^atfc] > q 
Dfc = <^ (A8) 
[ \ [(a'=at'=) - |(a'=)|2 - |Aa'=|2] if [a^at'=] < 

Eqn(A8) defines our fc*^ order squeezing parameter (Pfc), which has been used as the measure of squeezing in this 
paper. Clearly, 

Pfe < (A9) 
I 
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